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We work out the impact that the recently determined time-dependent component of the 
Pioneer Anomaly (PA), if interpreted as an additional exotic acceleration of gravitational 

'-fi ' origin with respect to the well known PA-likc constant one, may have on the orbital 

^ t motions of some planets of the solar system. By assuming that it points towards the 

If) ■ Sun, it turns out that both the semi-major axis a and the eccentricity e of the orbit of 

a test particle would experience secular variations. For Saturn and Uranus, for which 
modern data records cover at least one full orbital revolution, such predicted anomalies 

C^J , arc up to 2 — 3 orders of magnitude larger than the present-day accuracies in empirical 

determinations of their orbital parameters from the usual orbit determination procedures 
in which the PA was not modeled. Given the predicted huge sizes of such hypothetical 
signatures, it is unlikely that their absence from the presently available processed data 
can be attributable to an "absorption" for them in the estimated parameters caused by 
the fact that they were not explicitly modeled. The magnitude of a constant PA-type 
acceleration at 9.5 au cannot be larger than 9 X 10~ 15 m s — 2 according to the latest 
observational results for the perihelion precession of Saturn. 
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1. Introduction 

According to the latest analysis^ of ex tended data records of the Pioneer 10/11 
spacecraft, the small frequency driflP^ (blue-shift) observed analyzing the naviga- 
tional data of both the spacecraft, known as Pioneer Anomaly (PA), may present 
a further time-dependent component in addition to the well known constant one. 
Both linear and exponential models were proposed-^ for the PA; according to the 
authors of Ref. [H the exponential one is directly connected to non-gravitational 
effects^ since it takes into account the possible role of the on-board power genera- 
tors suffering a radioactive decay. 

In this letter we work out the orbital effects of such a new term in the hy- 
pothesis that the time-dependent PA component is due to some sort of long-range 
modification of the known laws of gravitation resulting in an additional anomalous 
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acceleration with respect to the nearly sunward constant one, having magnitude^ 

|Ap io | = (8.74 ± 1.33) x 1(T 10 m s~ 2 , (1) 

in terms of which the constant part of the PA has often been interpreted. Indeed, 
in this case it should act on the major bodies of the solar system as well, especially 
those whose orbits lie in the regions in which the PA manifested itself in its presently 
known form. In this respect, we will not consider the exponential model. Recent 
studie d I I I I I I I, partly preceding the one in Ref. [1] pointed towards a mundane 
explanation of a large part of the PA in terms of non-gravitational effects pertaining 
the spacecraft themselves. 

2. The orbital effects of the linear model 

Since the anomalous acceleration ia^ 

ipio ps -2 x KT 11 ms~ 2 yr _1 , (2) 

the time-dependent linear component of the postulated PA-type acceleration^ 

A = (t- t ) ipio (3) 

can be treated as a small perturbation of the dominant Newtonian monopole A^ 
over timescales of the order of an orbital period Pb for all the planets of the solar 
system. Table [T] explicitly shows this fact for Saturn, Uranus, Neptune and Pluto 
which move just in the spatial regions in which the PA perhaps started to appear 
(Saturn), or fully manifested itself (Uranus, Neptune, Pluto) in its presently known 
form. Thus, the Gauss equations for the variation of the osculating Keplerian orbital 

Table 1. Comparison between the magnitudes of the Newtonian 
monopole accelerations An of the outer planets of the solar system and 
their putative Pioneer-type accelerations A of eq. l(3j over timescales 
comparable to their orbital periods i\. Also the values of the plane- 
tary semi-major axes a, in au, and eccentricities e are displayed. The 
gravitational parameter of the Sun is GMq = 1.32712 X 10 20 m 3 s~ 2 . 



Planet 


a (au) 


e 


|A N | (ms" 2 ) 


Pb (yr) 


\A\ (ms" 2 ) 


Saturn 


9.582 


0.0565 


6 x 10- 6 


29.457 


6 x 10- 10 


Uranus 


19.201 


0.0457 


2 x 10 -s 


84.011 


2 x 10" 9 


Neptune 


30.047 


0.0113 


6 x lO" 6 


164.79 


3 x 10~ 9 


Pluto 


39.482 


0.2488 


4 x 10~ 6 


247.68 


5 x 10~ 9 



elements-^, which are valid for any kind of disturbing acceleration A, independently 
of its physical origin, can be safely used for working out the orbital effects of eq. 
(J3|). In particular, the Gauss equations for the semi-major axis a and eccentricity e 
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of the orbit of a test particle moving around a central body of mass M are 



# = ^{^^/ + ^T[c08/ + i(l-f)]} 



(4) 



they allow one to work out the rates of changes of a and e averaged over one orbital 
period -P D as 

f>-(i)f (?).* — 

In eq. (O (<Nf/dt) K are the right-hand-sides of eq. |4]) evaluated onto the unper- 
turbed Kcplcrian ellipse. In eq. (j4|) Ar, At are the radial and transverse components 
of a the generic disturbing acceleration A, p = a(l — e 2 ) is the semilatus rectum, 
n = y/GM/a? is the unperturbed Keplerian mean motion related to the orbital 
period by n = 27r/P D , G is the Newtonian constant of gravitation, and / is the true 
anomaly. Since the new data analysis^ does not rule out the line joining the Sun 
and the spacecrafts as a direction for the PA, we will assume that eq. ^ is entirely 
radial, so that Ar = A, At = 0. Using the eccentric anomaly E as a fast variable 
of integration turns out to be computationally more convenient. To this aim, useful 
relations are 

dt = ( 1 ~ e T c i osg ) dE, 

+ i f E — e sin E—En-\-e sin Eq\ / n \ 

* - *0 = { n ) , (6) 

sin/ 
As a result, a and e experience non- vanishing secular variations 



y'l-e 2 sin.E 
1 — e cos E 



I da\ _ A P q J e(2+e) 
\ dt I 



GM ' 

(7) 

/dex _ A P a 2 e (2+ e )(l-e 2 ) 

\dt/ ~ 2GM 

Notice that eq. ([7]) are exact in the sense that no approximations in e were assumed. 
Moreover, they do not depend on io- 

In order to make a meaningful comparison of eq. (J7|) with the latest empirical 
results from planetary orbit determinations, we recall that modern data records 
cover at least one full orbital revolution for all the planets with the exception of 
Neptune and Pluto. The author of Ref. fl~3l in producing the EPM2006 ephemerides, 
made a global fit of a complete suite of standard dynamical force models acting on 
the solar system's major bodies to more than 400,000 observations of various kinds 
ranging over At = 93 yr (1913 — 2006). Among the about 230 estimated parameters, 
there are the planetary orbital elements as well. According to Table 3 of Ref. [13j 
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the formal, statistical errors in a for Saturn and Uranus are 

c7 (EPM2006) =4)256m 



so that 



< PM2006) =40,294m, (8) 

6 



(EPM2006) . c _i 

a- — 4o m yr 

(EPM2006) , OQ _i ' 9 J 

cr^ = 433 m yr , 



can naively be inferred for their rates by simply dividing cq. © by At. The PA 
was not modeled in the EPM2006. It is important to remark that the figure for <t + 

quoted in eq. © was obtained without processing the radiotechnical observations 
of the Cassini spacecraft. According to eq. (0, the putative PA-induced secular 
changes of the semi-major axes of Saturn and Uranus are 



• (Pi°; 

% 



= 42,505 myr 1 

(10) 



4 Pio) ) =290, 581 my^ 1 



These are about 3 orders of magnitude larger than eq. ((9j) : even by re-scaling the 
formal uncertainties of cq. |5) by a factor of 10, the PA-type anomalous rates of 
cq. (|10[) would still be about 2 orders of magnitude too large to have escaped from 
a detection. Such conclusions are confirmed, and even enforced, by using the latest 
results published in Ref. [14] whose authors explicitly estimated secular changes of 
the semi-major axes of the first six planets with the EPM2010 ephemerides based 
on more than 635,000 observations of different types over At = 97 yr (1913-2010). 
The authors of Ref. HH whose goal was a different one, did not model the PA. They 
obtairfl 

.(EPM2010) -iq -1 /ii\ 

a^ = 13 m yr , (11) 



which is 4 orders of magnitude smaller than the predicted value of cq. (|10|) . 

Incidentally, we note that if eq. (j3]) acted on, say, the Earth, it would cause a 
variation in its semi-major axis as large as 

4 Pio) ) = 14.5 m yr- 1 (12) 



*A 35.9% correlation between dt and at was reported in Ref. 1141 
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corresponding to a shift of 

A<4 Pio) = 1,348.5 m (13) 

over Ai = 93 yr. The formal, statistical uncertainty in the terrestrial semi-major 
axis id^ 

CT (EPM2006) = 13 g m> (M) 

i.e. about 4 orders of magnitude smaller than eq. ([13")) . Even by re-scaling cq. (|14[) by 
a factor of 10, eq. (1T3"1) would still be 3 orders of magnitude too large. The EPM2010 



I; 



.(EPM2010) „ 1n _5 _i , 1t .% 

a y e '=2x10 m yr , (15) 



cphemcridcfo yicl 



i.e. 6 orders of magnitude smaller than the prediction of eq. (JT2J). 

In principle, it may be argued that cq. (|3|) was not included in the mathematical 
models fitted to the planetary data in the orbit determination process, so that its 
signature may have partly or totally been absorbed in the estimation of, say, the 
planetary state vectors: thus, the entire observational record should be re-processed 
by using ad-hoc modified dynamical force models explicitly including cq. (|3]) itself. 
However, this argument may have a validity especially when the magnitude of a 
putative anomalous effect of interest is close to the accuracy of the orbit determi- 
nation process: it is not the case here. Mor eover, as far as the constant PA term 
is concerned, independent analyse d I I ' in which it was explicitly modeled as 
a radial extra-acceleration acting on the outer planets substantially confi rmed its 
neat incompatibility with the observations, as suggesteqj in earlier studie d I I 

Concerning the possibility that the PA is present at Saturn, the authors of Ref. Q] 
yield an equivalent acceleration 

A^P = (4.58 ± 11.80) x lO" 10 m s~ 2 (16) 

for its constant term. It is known t hat a radial constant extra-acceleration causes a 
secular precession of the perihelio n ° > ^ ' ' ' given by 



^ Pi °) = Apio71 ^. (17) 

na 

The present-day accuracy in empirically constraining any anomalous precession Acb 
of the perihelion of Saturn from the INPOPlOa ephemerides, in which the PA was 
not modeled, id^ 

CAro*, = 0-65 milliarcsec cty -1 . (18) 



b A 0.6% correlation between dgj and agj was reported in Ref. 1141 

c Concerning the inner planets of the solar system, from an analysis of their determined orbital 
motions the authors of Ref. [3] pointed out that there was no room for a constant, sunward PA-likc 
exotic acceleration acting on them. 



April 5, 2012 1:4 WSPC/INSTRUCTION FILE Pioneer'2011'IJMPD 



6 L. Iorio 



Thus, eq. (HZ]) and cq. (HHJ) yield 



^Pio 



< -^L=a AA , = 9 x 1(T 15 m s" 2 , (19) 

Vl — e z <l 



which is up to 5 orders of magnitude smaller than eq. (|16|l . 



3. Conclusions 

If the time-dependent linear part of the PA were a genuine effect of gravitational ori- 
gin causing an extra-acceleration, it would affect the orbital motions of the planets 
with secular variations of their semi-major axis and eccentricity if it were oriented 
towards the Sun. The magnitude of such putative effects for Saturn and Uranus is 
up to 2 — 3 orders of magnitude larger than the current accuracy in determining 
their orbital motions from the observations without modeling the PA itself. A par- 
tial removal of a PA-likc signature may, in principle, have occurred in the standard 
parameter estimation procedure, so that the entire planetary data set should be 
re-processed with ad-hoc modified dynamical models explicitly including a hypo- 
thetical time-dependent PA-type acceleration as well. However, supported by the 
results of analogous tests for the constant form of the PA actually implemented 
by independent teams of astronomers with different purposely modified dynamical 
models, we feel that it is unlikely that such a potential removal of a PA-type sig- 
nature can really justify the absence of such a huge effect in the currently available 
processed planetary data. Concerning the existence of a constant PA-type acceler- 
ation at heliocentric distances of about 9.5 au as large as 1.6 x 10~ 9 m s -2 , latest 
results about the maximum allowed size of anomalies in the perihelion precession 
of Saturn yield a figure 5 orders of magnitude smaller. 
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